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ABSTRACT 

The Kelvin-Helmholtz cooling epoch, lasting tens of seconds after the birth of a neutron star in a successful 
core-collapse supernova, is accompanied by a neutrino-driven wind. For magnetar-strength 10'^ G) large 
scale surface magnetic fields, this outflow is magnetically-dominated during the entire cooling epoch. Because 
the strong magnetic field forces the wind to co-rotate with the protoneutron star, this outflow can significantly 
effect the neutron star's early angular momentum evolution, as in analogous models of stellar winds (e.g. Weber 
& Davis 1967). If the rotational energy is large in comparison with the supernova energy and the spindown 
timescale is short with respect to the time required for the supernova Shockwave to traverse the stellar progenitor, 
the energy extracted may modify the supernova shock dynamics significantly. This effect is capable of producing 
hyper-energetic supernovae and, in some cases, provides conditions favorable for gamma ray bursts. We estimate 
spindown timescales for magnetized, rotating protoneutron stars and construct steady-state models of neutrino- 
magnetocentrifugally driven winds. We find that if magnetars are born rapidly rotating, with initial spin periods 
(P) of ^ 1 millisecond, that of order 10^'- 10^^ erg of rotational energy can be extracted in ^ 10 seconds. If 
magnetars are born slowly rotating (P > 10 ms) they can spin down to periods of ^ 1 second on the Kelvin- 
Helmholtz timescale. 



Subject headings: stars: magnetic fields ■ 
gamma rays: bursts 

1. INTRODUCTION 



A successful supernova leaves behind a hot deleptonizing 
protoneutron star (PNS). This newly born, contracting, com- 
pact object radiates its gravitational binding energy in neutri- 
nos, which ablate matter from its surface via energy deposition. 
The primary heating mechanisms are the charged-current pro- 
cesses of electron and anti-electron neutrino absorption on free 
nucleons: i/^n pe~ and P^P ne'^- The neutrino energy (typ- 
ically 10 MeV) is deposited in the final-state leptons, which 
are hydrodynamically coupled to the nucleons. In this way, a 
thermal wind composed of nucleons and leptons is created. At 
the densities and temperatures surrounding the PNS at birth, 
photons are completely trapped and are advected with the wind 
(e.g. Duncan et al. 1986; Qian & Woosley 1996). 

This neutrino-driven outflow emerges into the supernova 
post-shock environment - the supernova shock itself moving 
away from the protoneutron star at a velocity of 10-30,000 km 
s"' (e.g. Burrows, Hayes, & Fryxell 1995). The duration of the 
PNS wind is set by the cooling or Kelvin-Helmholtz timescale 
(tkh ~ GM^/L^R, where is the total neutrino luminosity) for 
radiating away heat and lepton number and is typically of order 
tens of seconds (Burrows & Lattimer 1986; Pons et al. 1999). 
The kinetic luminosity of the wind is typically less than 10"^^ erg 
s"' and decreases sharply as the neutrino luminosity decreases 
(e.g. Qian & Woosley 1996) and the protoneutron star cools, 
and so the addition to the asymptotic supernova energetics is 
small on the scale of the canonical supernova energy, 10^' erg. 
Furthermore, the total mass ejected during the cooling epoch is 
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< lO""* M0, depending upon how the start of the wind phase is 
defined - minor in comparison with the total mass ejected in a 
typical core-collapse supernova. Despite this status as a mere 
perturbation to supernovae in both mass and energy, PNS winds 
have been the focus of considerable recent work. In particular, 
because of their intrinsic neutron-richness and association with 
supernovae, research on PNS winds has focused on the possibil- 
ity that they might be the as yet unidentified astrophysical site 
for production of the r-process nuclides (Woosley et al. 1994; 
Takahashi et al. 1994; Qian & Woosley 1996; Cardall & Fuller 
1997; Sumiyoshi et al. 2000; Otsuki et al. 2000; Wanajo et al. 
2001; Thompson et al. 2001). 

Recently, Thompson (2003 a,b) noted that in the absence of 
rotation magnetar-like (e.g. Kouveliatou et al. 1999; Duncan & 
Thompson 1992; Thompson & Duncan 1993) surface magnetic 
field strengths (Bq - 10^^ G) can dominate the thermal pressure 
and kinetic energy density of the wind in this very early phase 
of neutron star evolution. Here we extend this work by con- 
sidering in a simplified model the combined action of neutrino 
heating, rotation, and strong magnetic fields. 

Magnetic dipole radiation (e.g. Pacini 1967, 1968; Gunn & 
Ostriker 1969) and, to a lesser extent, gravitational wave ra- 
diation (e.g. Ostriker & Gunn 1969) are often considered as 
dominant spindown mechanisms, controlling the rate of rota- 
tional energy extraction and the angular momentum evolution 
of young neutron stars. Here we consider a mechanism well- 
known to the solar physics community: magnetic braking by 
the combined action of a persistent outflow and a strong mag- 
netic field (Weber & Davis 1967; Mestel 1968; Belcher & Mac- 
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gregor 1976; Hartmann & Macgregor 1982; Mestel & Spruit 

1987). For a wide range of initial rotation periods, we find that 
protoneutron stars may be spun down by the presence of this 
neutrino-magnetocentrifugally driven outflow in the first tens of 
seconds of the PNS's life. If these highly magnetic objects are 
born rapidly rotating, ~ 10^^ erg may be extracted on the cool- 
ing timescale, giving an energetic boost to the just-preceding 
supernova and providing favorable conditions for gamma ray 
bursts. 

Although our focus is on magnetars born with millisecond 
periods, we also assess the evolution of highly magnetic ob- 
jects bom slowly rotating, with spin periods similar to those of 
radio pulsars (tens to hundreds of milliseconds). Our conclu- 
sions for these objects are sensitively dependent upon magnetic 
field geometry. However, for monopole-like fields we come 
to the interesting conclusions that these objects are spun down 
considerably on the cooling timescale; a PNS with ~ 20 ms ini- 
tial period is transformed into a ^ 1 second rotator in tkh- 

1.1. This Paper 

In §2 we review the basics of stellar spindown via magne- 
tocentrifugal outflows and provide some motivation for pursu- 
ing this effect in PNS winds. Section 3 discusses some as- 
pects and ambiguities of the expected magnetic field topol- 
ogy. In §4 we provide simple analytic scalings for the Alfven 
point and spindown timescales over a wide range of initial 
PNS spin periods. Section 5 describes our models of neutrino- 
magnetocentrifugally-driven outflows and presents estimates of 
rotational energy extraction. In §6 we discuss applications to 
hypernovae and gamma ray bursts, making contact with Usov 
(1992) and Thompson (1994). We further summarize our find- 
ings, review the potential implications for supernova remnants 
and the angular momentum evolution of young neutron stars, 
and highlight open questions in need of further investigation. 
Throughout this paper we refer to models of non-rotating, non- 
magnetic PNS winds (e.g. the models of Duncan et al. 1986; 
Takahashi et al. 1994; Qian & Woosley 1996; Sumiyoshi et 
al. 2000; Otsuki et al. 2000; Wanajo et al. 2001; Thompson et 
al. 2001) with the acronym 'NRNM'. 

2. SPINDOWN 

Mass loss from the surface of any rotating star carries away 
angular momentum. If the star has a strong magnetic field, the 
matter lost in the stellar wind is forced into near corotation with 
the stellar surface out to ^ Ra, the Alfven point, where the 
magnetic energy density equals the kinetic energy density of 
the outflow. This effect provides for efficient angular momen- 
tum transport from the rotating star to the outflow (Schatzman 
1962). Angular momentum conservation implies that 
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surface Br B^, and the field is monopolar, Br oc r~^.^ Con- 
sideration of the azimuthal momentum equation together with 

Faraday's law yields 
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where p is the mass density. Although the magnetic field lines 
bend in the -0 direction between the stellar surface and /?a. 
such that BriRh) ~ B^{Rp,), eq. (2) states that the wind has an- 
gular momentum as if strict co-rotation is enforced between the 
magnetic field footpoints and 7? a- 

Employing eq. (2) in eq. (1) and taking / = (2/5)M/?^, where 
M is the PNS mass and is the PNS radius*, one finds that 

0/ = n,(M^/M,)'*^^/*^^ (3) 

Here we have assumed that R,^ and R^ are independent of time. 
The subscripts / and i denote 'final' and 'initial', respectively. 
In the absence of magnetic fields, R^ = Ry and for a total mass 
loss of even lO"'^ M©, changes by just 1% for a solar mass 
neutron star. For strong magnetic fields, however, we expect a 
large R^/Rv, and, correspondingly, large Ail. If strongly mag- 
netized PNSs are born rapidly rotating (Duncan & Thompson 
1992; Thompson & Duncan 1993) then a reservoir of rotational 
energy, 

-m^ ~ -I 

2 5 

large in comparison with 10^^ erg, may be extracted on the 

spindown timescale. In eq. (4), M1.4 = M/1.4 Mq, = 
R,j/\Q km, and Pi is the spin period in units of 1 millisecond 
(ms). The spindown timescale is defined here as the e-folding 
time for (combining eqs. 1 and 2); 



^Rot = ^/f^' ~ ^M/?202 ^ 2.2 X 10^2 erg M^^RI^^ F{\ (4) 



Tj = 



Q. 2M{R^ 



n 5 M \Ra 



(5) 



If Tj is short compared with the Kelvin-Helmholtz cooling 
timescale for the PNS, then the angular momentum of the PNS 
is significantly effected during the wind epoch. Also, because 
^Rot is potentially much larger than 10^' erg, if n is short com- 
pared with the time for the just-preceding supernova Shockwave 
to traverse the progenitor, then we expect interesting observa- 
tional consequences and modifications to the supernova explo- 
sion, resulting nucleosynthesis, and remnant dynamics. The 
last of these, that the supernova remnant may have larger in- 
ferred asymptotic kinetic energy is not unique to this spindown 
mechanism. Simple vacuum dipole spindown also energizes 
the supernova remnant on a relatively short timescale for highly 
magnetic, rapidly rotating neutron stars. The primary differ- 
ences here is that if a large fraction of the rotational energy can 
be tapped quickly, the nucleosynthesis of the actual explosion 
can be modified. This difference is of particular importance 
because, as we show in §4, the spindown timescale given in 
eq. (5) can be significantly shorter than that inferred from vac- 
uum dipole spindown. 



where / is the moment of inertia, is the angular velocity at 
the stellar surface, M is the wind mass loss rate, and C is the 
specific angular momentum carried by the wind. In the clas- 
sic model for solar spindown constructed by Weber & Davis 

(1967), the problem is treated in one spatial dimension at the 
equator with B = BrCr + B^e^ and v = Vr^r + v^e^. At the stellar 

' This field configuration is often referred to as 'split-monopole' . 

* The subscript u refers to the fact that in the PNS case we take the 'siuface' to correspond to the radius of neutrino decoupling, the neutrino-sphere 



3. FIELD TOPOLOGY 

To calculate the spindown of PNSs, we must estimate Ra- 
To do so, we appeal to its definition - the point at which the 
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magnetic energy density (B^/Sn) equals the radial wind kinetic 
energy density (pv^/2). For the purposes of simpUcity and gen- 
erality, we take 

Br = BoiR./rr. (6) 

The Weber-Davis model posits a monopole field structure with 
7] = 2. With this dependence of B on r, we may estimate and 
Tj. We may rightly ask, however, whether or not these estimates 
are appropriate considering the fact that the currents in the star 
should necessarily produce - to lowest order - a dipole field 
with 77 = 3. This question is germane because the predicted Ra 
is much smaller if one assumes a dipole field and, consequently, 
the inferred tj is much longer. The ambiguity is due to the in- 
trinsically multi-dimensional nature of magnetically-dominated 
outflows. Analytical (e.g. Mestel 1968; Mestel & Spruit 1987) 
and numerical models (e.g. Pneumann & Kopp 1971; Steinolf- 
son et al. 1982; Usmanov et al. 2000; Lionello et al. 2002; Ud- 
Doula & Owocki 2002) show that the interaction between the 
magnetic field and the flow is complex. The steady-state struc- 
ture, for modestly magnetically-dominated flows, is the classic 
helmet-streamer configuration (see Fig. 1 of Mestel & Spruit 
1987 for a depiction) with magnetic field lines emerging from 
low latitudes forming a 'dead zone', a region of closed mag- 
netic loops, and the wind at high latitudes opening the field to 
infinity. The resulting field is a mix of 'monopole' and 'dipole' . 

Because the kinetic energy density of the outflow dominates 
the magnetic field for r > Ra, and vice-versa for r < Ra, one 
might naively expect that inside Ra the field is roughly dipolar 
and outside Ra it is monopolar. This picture is the most conser- 
vative from the point of view of spindown and describes the ba- 
sics of the model of Mestel (1968). Rapid rotation complicates 
the issue. Mestel & Spruit (1987) show that the radial extent of 
the dipolar region, the dead zone (Rd), can be smaller than Ra 
(also impUed by Ud-Doula 2002). Assuming isothermal condi- 
tions, the radial extent of the dead zone can be calculated from 
the equations of magnetohydrostatic equilibrium; 
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in eq. (6). First, any net twist to the magnetic field lines yields 
77 < 3 and spindown is enhanced (Thompson et al. 2001). Sec- 
ond, the braking indices inferred from observations of pulsars, a 
system in which the pure dipole limit should strictly obtain, are 
inconsistent with 77 = 3. Third, the dead zone on the surface of a 
PNS may be periodically opened by neutrino heating (Thomp- 
son 2003), footpoint motion and shear (Thompson & Murray 
2001), or MHD instabilities. In short, 77 must be less than 3, 
but it is greater than 2. Both are limiting cases. The complica- 
tions of rapid rotation (eq. 7), relativity (we encounter relativis- 
tic outflows in §4), and neutrino heating make this ambiguity 
resolvable only with multi-dimensional numerical simulations 
of thermal magnetocentrifugal winds. 

In what follows, we present estimates for arbitrary 77 and scal- 
ings for 7; = 2 and 77 = 3. As we have just argued, our results 
for the monopole case hkely over-estimate the efficacy of spin- 
down (Poe, Friend, & Cassinelli 1989; Keppens & Goedbloed 
1999; Van der Hoist et al. 2002; but, see Pizzo et al. 1983). 
However, the results from the pure dipole (77 = 3) hmit must 
under-estimate spindown (Mestel & Spruit 1987). 



4. ESTIMATES 



At Ra, B^/Stt 
M/4TTr\, 



pv^/2. With eq. (6), and using p 



R 



2jj-2 . 



bIr^j^m- 
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(8) 



where va = Vr(/?A) is the Alfven speed. Let denote the 
asymptotic velocity attained by the matter in a purely NRNM 
outflow (that is, in a wind with f7 = and Bq = 0). Because the 
flow in NRNM wind models is driven by extremely inefficient 
neutrino heating, Vi, is typically much less than the PNS escape 
speed (see e.g. Thompson et al. 2001). For high neutrino lumi- 
nosity Vj, < 3 X 10^ cm s~^ If the matter is forced to co-rotate 
to ^ Ra, then v^(Ra) ^ Ra^- If Ra^ ^ i^V^ we expect that 
VriRA) ~ V0(/?a)- These arguments imply that to reasonable ap- 
proximation Va ^ Ra^- Therefore, from eq. (8), 



R^;>-'=BlRl^Ar'fi-K 



(9) 



where ct^ and denote the isothermal sound speed and mass 
density, respectively, in the dead zone, at R,, Equation (7) is a 
consequence of balancing magnetic tension with thermal pres- 
sure in an isothermal magnetohydrostatic atmosphere. Tak- 
ing fiducial PNS parameters as in §4 (/?^ = 10 km, P = 1 ms, 
Bo = 10'^ G, and M = 1.4 Mq), Rd is typically < 20 km, de- 
pending upon the choice for and ct„, much less than Ra 
(see eq. 10). In the range of radii Rd <r < Ra, the field fines 
are opened into a monopole configuration and Ra is estimated 
with 77 = 2 for the field structure. This effect yields spindown 
timescales that are shorter than those inferred from the more 
pessimistic case in which a dipole is assumed from the stellar 
surface all the way out to Ra- Therefore, eq. (7) implies that for 
rapid rotation (where 'rapid' is defined by the ratio /?^f2^/c| ) 
the net spindown might be better approximated by employing a 
pure monopole model, whereas for slow rotation a dipole field 
is probably more appropriate. 

In addition to this centrifugal effect, there are several other 
reasons to suspect that in estimating Ra we should take 7; < 3 

^ This is equivalent to assuming \hefast magnetic rotator (FMR) limit of Belcher iS 
* In the more general case (vi, ~ Ra^), one might take va ~ (vj +R\Cf-)^/'^ in eq. 



ScaUng for a rapidly rotating magnetar (P = 1 ms), 

RAiv = 2) ~ 43 km bI{^ Kio^Zl'^ P\'^ (10) 
RAir, = 3) ~ 24 km Bf^ R%1 ht'J^Pl" (11) 

Here and throughout this paper, Bo„ = Bo/10" G, M-m = 
M/IO"" Mq s"\ and Pi=P/l ms. Note that these scahngs ap- 
ply only when the condition Ra^ ^ is satisfied.^ 

There are several features of winds from rapidly rotating 
magnetically-dominated PNSs that are very different from their 
NRNM counterparts. Specifically, for the estimate given for Ra 
in eq. (9) to apply, Ra must be greater than /?sonic, the sonic 
point, where the radial velocity is equal to the local adiabatic 
sound speed (Cj). However, typical NRNM models yield sonic 
in the range of hundreds of kilometers, much larger than Ra as 
impUed by eqs. (10) and (1 1). This apparent problem is recon- 
ciled easily; in highly magnetic, rapidly rotating PNS winds, 
^sonic occurs at radii much smaller than those inferred from 
NRNM models. The boundary conditions at the PNS surface 

Macgregor (1976). 

(8). For comparison, Taam & Spruit (1989) employ va ~ [vi +(8/27)R^n^] V^. 
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- that the neutrinos are in thermal and chemical equilibrium 
with the matter and that the neutrino optical depth is ~ 2/3 
(see Thompson et al. 2001) - conspire to give a sound speed 
(Cs) at of ^ 3 X 10^ cm s~' . Therefore, for Riy = 10 km and 
spin period shorter than ~ 2 ms (fi > 3000 rad s"') Ri,Q ex- 
ceeds Cs(Ru)- This Umit, v^(R„) » Cs{Rv), is termed 'centrifu- 
gal' (Lamers & Cassinelli 1999). In this limit, the centrifugal 
force dominates wind driving and we may neglect both neutrino 
heating and the gas pressure in the critical wind equations near 
^sonic (see e.g. eq. 25; discussion of the general wind equations 
is deferred to §5). A simple expression results; 



(12) 



where Vg is the local escape velocity. The location of the sonic 
point is then 

= (GM/02)i/3 = {r^^rI/2)'I^ 



le.Skm Ml'] Pi^'\ 



(13) 

where /?sch is the Schwarzschild radius (~ 4.15 km for a 1.4 
Mq PNS). With this estimate we find that the sonic point oc- 
curs much closer to the PNS and that Ra > ^sonic for eqs. (10) 
and (11). 

There is another feature of centrifugal winds that differs sig- 
nificantly from NRNM outflows. As R,^n becomes greater than 
CsiRu) the matter density scale height increases exponentially 
with fp- due to centrifugal support in the near-hydrostatic re- 
gion (vr <C Cs). For this reason, the derived mass loss rate also 
increases exponentially with fi^. Note that this is only true if 
the magnetic field supports the flow from the surface of the star 
to the sonic point. That is, we only expect large increases in 
M for Ra > /?sonic- Examining eq. (5) we see that with M in- 
creasing exponentially with 17^, tj may become very short for 
large fl. However, as fl is increased and the mass loading on 
the field lines increases, Bo must be large enough to guarantee 
Ra > ^sonic (eq. 9). In NRNM models, characteristic mass loss 
rates are less than 10"^ Mq s"^ The exponential effect on M 
in the centrifugal limit justifies our use of a much higher mass 
loss rate in eqs. (10) and (11), when f = 1 ms (see Fig. 2 and 
§5). 

In an effort to make contact with the pulsar and solar wind 
communities, it is worth noting that taking eq. (9) and multi- 
plying both sides by fi^ for = 2, one obtains 

vi = vii^— f^ = -^, (14) 
M Fm 

where Fb is the magnetic flux, Fm is the matter flux, and vm 
is the Michel velocity (Michel 1969; Belcher & Macgregor 
1976). As one would expect from the Bernoulli integral, it is the 
Michel velocity that, to constants of order unity, is the asymp- 
totic wind velocity. In a true fast (vi, ^ Ra^), magnetically 
dominated wind the Alfven velocity is 2/3 of the Michel ve- 
locity, and Vm obtains only at the fast magnetosonic point (e.g. 
Belcher & MacGregor 1976). In fact, for a strictly monopolar 
cold (cj = everywhere) magnetocentrifugal wind the fast point 
is formally at infinity (Michel 1969). At the level of accuracy 
aspired to here, these details are unimportant. 

Spindown timescales (tj = f2/f2; eq. 5) can be computed 
from the estimates of Ra in eqs. (10) and (11). Together with 
eq. (9) we find that 

= (2/5)MA^3-2,)/(2„-l)^--2/(2,-l)g^/(2^-l) j^2/(2^-l)_ (^5) 



The superscript 'NR' is added to emphasize that when the flow 
is non-relativistic, rj depends explicitly on M as in eq. (5). For 
relativistic flows it does not (see eq. 22). Taking M=\A M© 
and scaling, 

^\r^ = 2) ^ 30 s M,.,M-_l/'Riy'B^fpf' (16) 
TfHv = 3) ^ 96 s M^.,M:l/'R;fB^fpf' (17) 

The important timescale to compare with rf"* is the Kelvin- 
Hehnholtz cooling timescale, tkh- This comparison is relevant 
because M drops steeply with the PNS neutrino luminosity. If 
Tj^ ':$> Tkh no spindown occurs during the wind epoch. The 
PNS cooling calculations of Pons et al. (1999) show that 
drops by a factor of ten, from ^ 10^^ erg s"' to ^ 10^' erg s"\ 
in ~ 30 seconds. This result depends on the mass of the PNS 
and the high-density nuclear equation of state employed. These 
calculations do not include the effects of rapid rotation or high 
magnetic fields. Both may significantly modify tkh- Rapid 
rotation leads to lower core temperatures and, thereby, lower 
average neutrino luminosity. This could increase tkh signifi- 
cantly (see Thompson, Quataert, & Burrows, in prep; Yuan & 
Heyl 2003; Villain et al. 2003). Rapid rotation might also make 
Ri, contract more slowly, affecting both Bq and 51 at a given 
time because both are proportional to R~^, via flux and angular 
momentum conservation, respectively. High Bq may also affect 
Tkh by modifying neutrino opacities (e.g. Lai & Qian 1998; Ar- 
ras & Lai 1999). Multi-dimensional effects such as convection 
may also be important (Keil, Janka, & Miiller 1996). With- 
out detailed models of magnetar cooling, we take the results of 
Pons et al. (1999) as representative and define tkh ~ 30 sec- 
onds. 

As argued in §3 the spindown timescales relative to tkh are 
sensitive both in absolute value and scaling to the parameter rj. 
Even allowing for our ignorance about the magnetic field topol- 
ogy, eqs. (16) and (17) imply that at best t^^ ~ tkh for the PNS 
parameters chosen. Note that for higher M, both timescales de- 
crease somewhat, despite the fact that Ra decreases as M in- 
creases. Interestingly, for larger initial spin periods, the spin- 
down timescales decrease. Hence, if magnetars are born with 
P 30 ms - consistent with the spin periods of most pulsars 
today - then t^^ ^ 3 s and ^ 25 s for 77 = 2 and 77 = 3, respec- 
tively. As we will see in §5 mass loss rates as high as 10~^ M© 
s"' are probably unrealistic at such low fi. Even so, spindown 
of slowly rotating magnetars during the cooUng epoch can be 
quite efficient. For higher Bq, both spindown timescales drop 
rapidly. Increasing the surface magnetic field strength to lO'^ G 
in eq. (17), rj^^irj = 3) ~ 15 s, which may be considerably less 
than Tkh- Although such a large is not out of the question 
in this very early phase of PNS evolution, it is well beyond the 
surface field strengths inferred from magnetars (Kouveliotou et 
al. 1999; Thompson & Duncan 1993). 

Consideration of such high field strengths brings up an im- 
portant physical constraint on Ra, particularly germane for 
PNSs rotating rapidly at birth: Ra cannot be made arbitrarily 
large by increasing Bq for a given Q because would eventu- 
ally exceed the speed of light. Thus, the Alfven point, deter- 
mined a posteriori for a given il, must be less than the light 
cylinder radius, Rl = c/fl, although for asymptotically large ra- 
tios of magnetic flux to matter flux Ra asymptotes to This 
requirement sets a critical Bq. If we set Ra = Rl = c/ft in eq. (9) 
for a given fi, R^, and M, the Bq derived is that required to en- 
force co-rotation to approximately R^. This critical magnetic 
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field, for arbitrary r], is easily derived in terms of the basic wind 
and PNS parameters; 



Br 



(18) 



For example, for a neutron star with spin period of 1 ms, 
/?L 47.7 km and for M = 10"^ M© s~^ as in eqs. (16) and 
(17), Bcrit ^ 1.2 X 10^^ G and 5.7 x lO" G for ry = 2 and r/ = 3, 
respectively. 

If Ra is sufficiently less than the asymptotic velocity of 
the wind material is not relativistic. In non-relativistic out- 
flows, consideration of the Bernoulli integral in the Weber- 
Davis model shows that the energy carried by the magnetic field 
is a factor of two greater than that of the matter (e.g. Lamers & 
CassinelU 1999) and the rotational energy loss rate is 

£,ot = mn = -Mcn = -B^^E^ji'^M^-^'^n'^-^^, (i9) 

where C = 1 /(2t?- 1). If Bo exceeds Bcrft, approaches /?l and 
va approaches c. The flow becomes relativistic and the field 
may carry more than two times the energy of the matter. At 
the ratio of magnetic flux to matter flux is given by^ 

" = BlRl^n'''>-^c'-^'>i\rK (20) 



and the relativistic energy loss rate is 

£,ot = -FMc^ = -BlRl^^f^-V-^\ (21) 

Note that Etot is independent of M, but depends on the magnetic 
field structure. For the case r] = 3, the pure magnetic vacuum 
dipole scaling for fitot is obtained.^" For = 2, this expression 
for Etot implies a much larger energy loss rate than for rj = 3 - 
the ratio being c^/Rln^, a factor of ~ 23 for a 10 km PNS wifli 
a 1 ms spin period. As M drops during the cooling epoch, for 
constant Bo, the wind transitions from a non-relativistic (eq. 19) 
to a relativistic outflow (eq. 21). Only in the case where f2 is 
very small so that /?l is very large are there instances when 
the wind is always non-relativistic during tkh (see §5.3). Tak- 
ing eq. (20), the spindown timescale (rj = Q/Q) for relativistic 
winds, in analogy with eq. (15), is 

(22) 



4-2jj^2))-3 



rf = (2/5)M/?2-2''B„2f^4- 
Taking the same scalings as in eqs. (16) and (17) 
Tf(r? = 2): 



34sM,.4/?;iBn' 



Tf(r, = 3) ~ 760 sMi.4«^i 



(23) 
(24) 



of eq. (24) when the mass flux is ~ 10"^ Mq s"' under-predicts 
the rate of angular momentum loss by a factor of 8." Note 
that this discrepancy is much larger for slowly rotating magne- 
tars; Tf-(ri = 3)/T^^(ri = 3) oc f '2/5_ jf magnetars are born with 
P = 10 ms, Tj^/tj^** ~ 2000 and one may under-predict spin- 
down by several orders of magnitude. Of course, the difference 
is further magnified if, when the mass loss rate is high, the field 
is monopolar to some extent, as argued in §3. In this case we 
compare eq. (24) (rj = 3) with eq. (16) (rj = 2) and find that the 
vacuum dipole limit overestimates the spindown timescale by 
a factor of ^ 25 for the fiducial 1 ms rotator with Bq = 10'^ G. 
Noting that (?? = 3)/rf "^(r? = 2) cx f we see that this ratio 
is ~ 10"^ for a slow rotator with P = 10 ms. 



5. NEUTRINO-MAGNETOCENTRIFUGAL WINDS ON THE 
CHEAP 

Because the spindown timescales and the amount of rota- 
tional energy extractable depend so crucially on M, we solve 
here the steady-state wind equations including neutrino heat- 
ing and rotation explicitly. Magnetic fields are included impUc- 
itly by enforcing = rfi everywhere. Thus, the wind mass 
elements are forced to corotate with the PNS surface as if on 
wires emanating radially from the equator of the PNS like the 
spokes on a bike wheel. This approximates the effect of a strong 
magnetic field or the non-relativistic Weber-Davis problem. Of 
course, in our solution the mass elements eventually have ve- 
locity greater than c. For a given O, the radius at which this 
occurs is the light cylinder For this reason our solutions are 
valid (if approximately) only for r < Importantly, since it is 
the effect of the centrifugal force on M that we are most inter- 
ested in, and M will be most effected by the centrifugal force in 
the subsonic region, the sonic point (see eq. 13) is always well 
within /?L where the corotation assumption is valid. 

It would be preferable to solve the one-dimensional relativis- 
tic Weber-Davis problem (Michel 1969; Goldreich & Juhan 
1970) with an arbitrary energy deposition function. We save 
this for a future paper. Much more important, we feel, is our 
assumption that the problem may be treated in one spatial di- 
mension (as in Weber & Davis 1 967) in the presence of rapid ro- 
tation and strong magnetic fields (see §3). This approximation, 
can only be removed with a code capable of two-dimensional 
axisymmetric relativistic MHD for /? <C 1 . This is a consider- 
able effort beyond the scope of the current work. 



5.1. Equations 

Enforcing continuity, conservation of momentum, and con- 
servation of energy in Newtonian gravity, we have that 



The superscript 'R' is meant to reinforce the difference between 
relativistic and non-relativistic wind spindown - the latter de- 
pending explicitly on M. Comparing eqs. (24) and (17) we 
see that rf^** is significantly shorter than rf. The difference is 
one of applicability. Equation (17) applies with large M, when 
Ra < Rl and F < 1 . Equation (24) applies only when the mass 
flux from the PNS subsides to such an extent that F > 1 . Tak- 
ing Bo = Bcrit in eq. (17), t^^ = rf-. Therefore, application of 
vacuum dipole spindown is appropriate only when M is suffi- 
cienfly small. We see that a naive and inappropriate application 

' The ratio F, as defined in eq. (20), is often denoted by 'a' in the pulsar literature. 

Note that in the standard theory (see e.g. Shapiro & Teukolsky 1983 chapter 10 for a review) Em in the vacuum dipole limit includes a factor sin^ a/6, where a is 
the angle between the spin axis of the neutron star and its magnetic axis. 

" Although both considered magnetars in the context of soft gamma repeaters and anomalous X-ray pulsars, Thompson & Blaes (1998) and Harding et al. (1999) 
showed that enhanced spindown is obtained in models including relativistic particle winds. 

We have developed a one-dimensional non-relativistic time-dependent Eulerian magnetohydrodynamics code for solving the isothermal Weber-Davis problem. 
Preliminary comparisons between models in which strict corotation is enforced (v^ = rQ) everywhere and the actual Weber-Davis solution show that Sa is typi- 
cally ~ 1 .5 times larger than that expected from the models presented here. This implies that we underestimate the magnitude of the spindown for a given Bp and 
overestimate spindown timescales r, by a factor of ~ 2. 



Vr dr ^ ' r\A 2 



+ q 



D 

CvTVr 



(25) 
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,26) 

par r \ 4 2 1 CyiVr 

(27) 

where v, = (2GM/r)i/^ z) = (T/p)dP/dT\p, ct and c, are the 
isothermal and adiabatic sound speeds, respectively, Cy is the 
specific heat at constant volume, = ril, Vy is the radial ve- 
locity, and q is the specific neutrino heating rate. For a given 
f2, we solve the above equations for the flow between the sur- 
face of the PNS and the sonic point as a two-point boundary 
value problem, using a relaxation algorithm, on an adaptive ra- 
dial mesh (see Thompson et al. 2001). We integrate from the 
sonic point to the light cyUnder using a simple Runga-Kutta al- 
gorithm. The equation of state (EOS) assumes ideal nucleons. 
It includes a general electron-positron EOS, photons, and ne- 
glects the formation of alpha particles. Magnetic effects on the 
EOS, particularly important for electron-positron component at 
low temperatures, have not been included. Once a solution is 
obtained, we can see what surface magnetic field strength is 
required to corotate out to any radius. Once the field is cho- 
sen, the solution is vahd only inside this a posJenon-posited 
Alfven point. Outside this radius, our solutions continue to ac- 
celerate as a result of the centrifugal force. In the full Weber- 
Davis problem, begins to decrease at Rp, and v,. asymptotes, 
reaching the Michel velocity at the fast magnetosonic point (see 
Belcher & Macgregor 1976). 

5.2. NMC Wind Results 

Figure 1 shows the thermal pressure (solid line) and kinetic 
energy density (pv^/2 long dashed line and pvl /2 short dashed 

line) for a PNS with Lp^ = 4 x 10^' erg s"', M = 1.4 Mq, and 
Ri, = 10 km. We use L^^ to label wind models; the total neu- 
trino luminosity for a given model is L'°' = Li,^ +Li^^ +4Li,^, 
with L[y^_ = L[yJ1.3 and L^^ = Lp^/1.4. The subscript V^' 
stands for i^^, i/^, Pt, and Dt- For L^^ = 4 x 10^^ erg s~\ 
1}°'- ~ 1.85 X 10^^ erg s~K The PNS period is taken to be 1 ms 
(fi = 6283 rad s"'). Overlayed on Fig. 1 is the magnetic energy 
density (dotted lines, labeled B^/Stt) for rj = 2 and rj = 3 such 
that B^/Sn ^ pc^/2 at Ri^. The point at which the dotted curves 
intersect the short dashed Unes corresponds to the Alfven ra- 
dius. The surface magnetic field strengths are labeled Bq and 
correspond to Bent in eq- (18) for q = 2 and 77 = 3. For these 
field strengths, then, the Alfven point corresponds to ~ For 
lower Bo, Ra decreases. Figure 1 shows that we can expect a 
1.6 X 10^^ G surface magnetic field with monopolar field topol- 
ogy to enforce corotation out to ^ /?l for mass loss rates in the 
range of 10"-' M0 s"' (the model shown has M ~ 1 .6 x 10"^ M0 
s"^). For = 3, the required surface field strength is ~7. 5x10'^ 
G. As L^°^ increases M increases and larger Bq is required to 
force corotation out to R]^. As L'°' decreases, the converse is 
true. For example, for Lp^ = 0.5 x 10^^ erg s~\ the required sur- 
face field strengths are L4 x 10^^ G and 2.9 x lO"* G for 7/ = 3 
and 77 = 2, respectively. That B^nt decreases as L'°' decreases is 
expected from eq. (18) and the fact that the mass loading of the 
field lines, M, decreases steeply with L!^^. 

The intersection of P with pvf/2 (at r 17 km) marks the 
radius of the sonic point and is very well described by eq. (13) 
for short PNS spin periods, regardless of L'°'. Note that for 



the Bo used to plot B'^/Stt in Fig. 1, for both r/ = 2 and = 3, 
that deep in the exponential atmosphere of the PNS, both P 
and pv^/2 become greater than the magnetic energy density. 
This is potentially important since we have assumed in con- 
structing these models that corotation obtains (v^ = r^l) every- 
where. This assumption is most important between R^, and 
^sonic, because centrifugal support in this region can increase 
M significantly. The surface of the PNS calculations, 7?^, is 
defined as the neutrinosphere (optical depth ^ 2/3). This con- 
dition determines p{Rv), which, for most calculations is ^ IQ^^ 
g cm"-'. The azimuthal kinetic energy density at Ri, is then de- 
termined by setting il(R,y). Although in Fig. 1 and §4 we have 
assumed that R^, corresponds to the radius at which the mag- 
netic field is anchored (B(r) = BoiRv/r)'^), this need not hold 
necessarily. In addition, very close to the PNS, where the radial 
flow is very subsonic, we expect the field to be a complex of 
higher-order multi-poles as shear energy within the PNS con- 
vective region emerges as magnetic flux into the magnetically- 
dominated atmosphere via the Parker instabiUty (Thompson & 
Murray 2001). The degree of magnetic support within this at- 
mosphere, particularly if differentially rotating, will depend on 
the dynamics of the underlying convection and the efficacy of 
the dynamo mechanism, issues that require a multi-dimensional 
model of these MHD winds. 

The importance of the combined action of rapid rotation and 

strong magnetic fields is seen most clearly in Figure 2. Here, 
we present the mass loss rate M for models with Lp^ = 8 x 10^^ 
erg s~' and 0.5 x 10^' erg s"^ as a function of fi. This range of 
neutrino luminosity corresponds to roughly 20 seconds of the 
cooUng epoch in the models of Pons et al. (1999). The soUd 
lines assume strict corotation with the PNS surface, = rfl. 
The dashed lines are presented for comparison assuming no 
corotation and angular momentum conservation; in these mod- 
els = Ri,n{Ri,/r). All mass loss rates quoted are 'spherical' 
and assume M = 4TTr^pVr. For il > 3000 rad s"' in the models 
assuming corotation, significant deviations from the = limit 
are obtained. This results from the fact that at this for Ri,= lO 
km, Ri,fl becomes greater than CgiRv). As alluded to previ- 
ously (§4), this changes the mass flux in a dramatic way; for 
n = 6283 rad s"' and Lp^ = 8 x 10^' erg s'^ M ~ 8 x 10^^ Mq 
s"', more than 30 times larger than in the non-rotating model 
(M ~ 2.6 X 10-4 Mq s"^). For comparison, wifli n = 9000 rad 
s-\ M ~ 0.36 M0 s-\ For Q > 3000 rad s"', the functional 
dependence of M can be approximated roughly with 

MiflR^ > CsiR^)) oc exp [Tl^Rl/c^^iR^)] (28) 

as expected from simple considerations of the subsonic expo- 
nential atmosphere of the PNS. For fl ~ 3000 rad s"', M de- 
pends on the location of the sonic point itself (see e.g. Belcher 
& Macgregor 1976; Hartmann & Macgregor 1982; Lamers 
& Cassinelli 1999). For rapidly rotating PNSs the spindown 
timescale may be very short as a result of the large increase 
in M expected when R^Q Cj (eqs. 16 and 17). As an ex- 
ample, taking Lp^ = 8 x 10^' erg s~', with P = I ms, r] = 2, 
and Bo - B,a - 3.6 x lO'^ G, we have from eq. (16) that 
rf ~ rf"* ~ 3.1 seconds. These numbers, together with the 
scaling relations for t^^ and in §4, imply that a significant 
amount of rotational energy may be extracted from the PNS 
during the neutrino cooling epoch. 
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5.3. Rotational Energy Extraction & Spindown 

The spindown timescale depends crucially on r/. This is seen 
in eqs. (15) and (22). In addition to the arguments of §3, there 
are a few things that can be argued with some certainty about 
the field structure. First, if B^/Stt <C P at the surface, then no 
closed magnetic field lines will persist and Br oc r"^, despite the 
fact that the currents in the star should necessarily produce a 
dipole field. Although uninteresting from the perspective of 
spindown because no Ra exists outside 7?^, certainly in this 
Umit we should have 77 ~ 2. Second, and at the other extreme, 
in a totally magnetically dominated system with F very large, 
we expect the field to be dipolar with closed field lines within 
Rl (although, see §3). In this hmit, if we accept the magnetic 
dipole spindown model for magnetars (and neutron stars gener- 
ally), then we must have 77 ~ 3. These simple arguments sug- 
gest that T] itself is a function of the ratio of the magnetic flux 
to the matter flux. Following Ud-Doula & Owocki (2002), one 
may construct the dimensionless number 

C = bIrI/Mvoo (29) 
and write r] as a function of ^. However, with the complica- 
tions of rapid rotation and relativity, we have decided to simply 
bracket the possible solution space, computing solutions with 
r] = 2 and r] = 3 separately as in §4. 

Monopole Spindown: We solve the differential equation 

/r217 = Etoi for n(t) with rj = 2. If Ra < Rl (non-relativistic 
wind), we use Etot = -MR\Q^, with Ra from eq. (9) with 
rj = 2. If r > 1 (relativistic wind), we take £tot = FMc^, with 
r = Bf^R'ln^M'^c'^ . In order to obtain M, we employ the re- 
sults from our NMC models. We assume that these steady- 
state wind solutions, for different PNS characteristics, may be 
concatenated to form a time series. The total neutrino heating 
rate is computed self-consistently, assuming a PNS neutrino lu- 
minosity, including contributions from z/^n ^ pe~, v^p ^ ne^, 
vv ^ e^e~, and inelastic scattering of neutrinos with charged 
leptons and nucleons (Thompson et al. 2001). For the models 
presented here, as reflected by the requirement that Ra^ is al- 
ways greater than (see §4), the total energy deposition rate 
via neutrino heating and cooling is more than an order of mag- 
nitude less Mv^. This implies that all wind models are pri- 
mary driven magnetocentrifugally. In NRNM models of PNS 
winds, if oc lU^, where is the average neutrino energy, 

M oc lI^^ (Qian & Woosley 1996).i3 With this dependence in 
hand from our wind models we need only posit L,y{t). Tak- 
ing L[y(t < 40 s) OC f"", with a in the range ^ 1 as appropri- 
ate over much of the cooling epoch (see Pons et al. 1999), and 
L^it > 40 s) oc with r = 1 second, we can compute flit). 

As an example of slow magnetic rotator spindown we take an 
initial spin period of 20 ms and Bo= 1.5 x 10'^ G. Integrating 
from ti = l s to tf = 100 s, we find that the PNS stops rotating 
in ^ 50 seconds. In this model a total of 1.75 x 10"^ Mq is 
ejected and all of the rotational energy of the PNS is extracted. 
For this choice of parameters, F < 1 throughout the evolution. 
For longer spin periods we find similar evolution. The PNS 
simply stops rotating in tens of seconds. Hence, if magne- 
tars are born slowly rotating (with an initial spin period simi- 
lar to normal neutron stars) we find that monopolar neutrino- 
magnetocentrifugal winds can efficiently extract the rotational 
energy of the young PNS on the Kelvin-Helmholtz timescale. 

This dependence of the mass loss rate on luminosity is only approximate. For 
and the exponent decreases. 



As an example of an evolutionary sequence in which we start 
with a rapidly rotating PNS, we take the initial spin period to 
be 1 ms and the same temporal evolution of the neutrino lumi- 
nosity, folding in the exponential dependence of M on il for 
n > 3000 rad s'^ (as in Fig. 2). Results for 5o = 3 x 10^^ G 
(A), 1.5 X 10^5 G (B), 10'^ G (C), and 7 x lO^'' G (D) are pre- 
sented in Figure 3. The posited time evolution of is shown 
in the lower right panel. The evolution of the mass loss rate 
follows from the models constructed in §5 and is shown in the 
upper right panel. Computed evolution of Pit) is in the upper 
left panel. For clarity of presentation, only the first ~ 20 sec- 
onds of Pit) for trajectory (A) are shown. The lower left panel 
shows the total rotational energy extracted from the PNS as a 
function of time (solid lines, 10^' erg) and the rotational energy 
loss rate (dashed lines, 10^' erg s"'). The lower middle panel 
shows r as a function of t for models (A)-(D), and a model 
with Bo = 1-3 X 10^^ G (heavy dotted line). Finally, the upper 
middle panel shows the time evolution of Rl (solid lines) and 
Ra with 7] = 2 (dashed lines). 

In model (A) the flow begins with Ra ^ Rl- For model (D) 
several seconds pass before the flow becomes relativistic and T 
becomes greater than 1 . For all models shown here, over the 
100 seconds of evolution computed, a total of > 2 x 10^^ erg 
in rotational energy is lost by the PNS. Because of our assump- 
tion of rj = 2, the PNS spin period continues to e-fold on the 
timescale given by eq. (23). The total mass lost by the PNS de- 
creases as B() increases; for model (A) 3.3 x 10"^^ Mq is ejected, 
whereas for model (D) 5.3 x 10"^ Mq is emitted. Because fl 
drops so rapidly in model (A), and M depends exponentially on 
n for P < 2 ms, this correlation is expected. In models (B), (C), 
and (D), more than 10^^ erg is emitted with F > 100. In fact, 
at the end of these spindown calculations, with M plummeting 
with the exponential cutoff in Lj^, F can easily reach 1000. The 
exact evolution depends sensitively on Bq and the time evolu- 
tion of M, as can be seen from the model with Bq = 1-3 x lO'^ 
G (dotted line) in the lower middle panel. The complex time 
evolution of F comes both from the exponential cut-off in M 
(imposed at ? = 40 s) and the fact that F can decrease even as 
the mass loss abates if Rl increases rapidly. That is, for con- 
stant B{), r can decrease if the spindown timescale is shorter 
than the timescale for M to decrease (see eq. 20). The non- 
relativistic epoch of spindown (F < 1) lasts a mere ~ 1 - 2 s 
in models (B)-(D), and accounts for ~ 2-3 x 10^^ erg in rota- 
tional energy extracted. Although this timescale is short on the 
scale of tkh, it is very long on the scale of the PNS dynamical 
time and the spin period. Model (A) exhibits begins with near- 
relativistic flow and only after r ~ 20 s does the flow become 
non-relativistic, with F peaking at ^ 100 at f ^ 6 seconds. We 
conclude that spindown in a Weber/Davis-type magnetic field 
geometry, taking rj = 2, can be significant during the cooling 
epoch. Even if a monopole-Uke field structure only obtains 
during the F < 1 non-relativistic epoch, an amount of energy 
comparable in magnitude with the supernova explosion itself 
(> 10^^ erg) is extracted and emerges into the post-explosion 
supernova ejecta. 

Dipole Spindown: We solve for flit) in the same way as in 
the monopole case, but here take 77 = 3 in calculating Ra and, for 
F > 1 use Eiot = BqP^JI'^c"^. For a slow rotator with an initial 

large Q (above 7000 rad s"') this power law dependence of M on L^, is modified 
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period of 20 ms and Bq = 10^^ G, the final spin period is ~ 21 
ms. That is, a lO'^ G dipole field does not give significant spin- 
down. For this slow initial rotation period, only a dipole field 
greater than lO'^ G affects the evolution significantly before 
r > 1 and the simple vacuum magnetic dipole hmit obtains. As 
an extreme case, taking Bq = 5 x 10^* G, the final spin period is 
^ 40 ms, more than 70% of the rotational energy is extracted 
(a mere ~ 4 x lO"*** erg), and T is greater than unity for / > 35 
seconds. Figure 4 shows the evolution of Pit), ^a, T, and 
E and E. Virtually all of the braking occurs in these first 35 
seconds. As the wind transitions from non-relativistic to rela- 
tivistic, spindown ceases on the timescale tkh and the dipole 
limit obtains. 

Figure 5 shows the time evolution computed for PNSs with 
initial spin period of 1 ms and employing 77 = 3. For compar- 
ison, the same evolution is shown for three different surface 
magnetic field strengths, Bo = lO^'^G, 5 x lO^^G, and lO'^G, 
labeled 'A', 'B', and 'C, respectively. Panels, linestyles, and 
units are the same as in Fig. 4. For Bo = lO'^ G, the final spin 
period is 1.15 ms. Although the total change in spin period 
is small, approximately 10^' erg is extracted in the first sec- 
ond of the spindown calculation. A total of 2 x 10^' erg 
is ejected during the non-relativistic phase of wind evolution 
(f < 9 s), with Ra < Rl and F < 1 . Of course, for the remain- 
ing 90 seconds of the spindown calculation, the simple vacuum 
dipole limit obtains. For higher field strengths the transition 
from non-relativistic (F < 1) to relativistic (F > 1) occurs at 
earlier times. For the highest field strength shown here (A), the 
magnetic energy density dominates the kinetic energy density 
at /?L throughout the evolution. Most important are our results 
from the (C) calculation, which show that even taking the more 
pessimistic dipole field (rj = 3 instead of r/ = 2), we find that 
an amount of energy comparable to the asymptotic supernova 
energy is naturally injected into the post-shock supernova envi- 
ronment on a short timescale with respect to the amount of time 
needed for the preceding supernova Shockwave to traverse the 
progenitor (see §6.1). 

Note that by decreasing the initial spin period, while keeping 
Bo = 10^^ G, the spindown timescale drops sharply as a result of 
the exponential dependence of M on f2 (see eqs. 9 and 28). For 
example, taking an initial period of 0.8 ms and Bo = 4 x lO'^ 
G, one finds M ~ 0.074 M© s"' and Ra ^ 17 km. This im- 
plies a non-relativistic (F < 1; see eq. 17) spindown timescale 



I J 2.6 seconds. For comparison, the vacuum dipole spin- 
down approximation (inappropriate with such large M) yields 
~ 31 seconds.'"* In the example given, with initial period 
0.8 ms and Bq = 4 x lO'^ G only three seconds transpire be- 
tween the start of the calculation and the moment the vacuum 
dipole limit obtains. However, in those few seconds, ~ 10^^ erg 
is extracted - a total of 0.019 M© with velocity ^ 1.5-3 x 10 



6. DISCUSSION 

6.1. Hyper-Energetic Supernovae 

For the purposes of this paper we define the word 'hyper- 
nova' as any supernova with an asymptotic observed total en- 
ergy of greater than 10^' erg. We propose that the magneti- 
cally dominated, centrifugally slung winds from rapidly rotat- 
ing PNSs can explain hypemovae. 

In this scenario, core-collapse is followed by shock stall and 
then explosion. The explosion may be caused by a combina- 
tion of neutrino heating and convection (e.g. Herant et al. 1994; 
Burrows, Hayes, & Fryxell 1995; Fryer & Warren 2002), by 
neutrino heating plus rotation and dissipative effects (Thomp- 
son, Quataert, & Burrows in prep) or perhaps via MHD-driven 
jets (Akiyama et al. 2003). In any case, the supernova shock 
propagates outward at ~ 10, 000 km s"' with energy ~ 10^' erg. 
The wind phase conraiences as the PNS contracts and cools. In 
the first few seconds, the kinetic energy density and thermal 
pressure of the wind are large and the wind is not magnetically 
dominated. The wind hydrodynamical power in this phase is 
inefficiently supphed by neutrino heating as in NRNM wind 
models. As the radius of the PNS decreases, the average mag- 
netic field grows because of flux conservation (Bo oc R~^) and an 
efficient dynamo as in Duncan & Thompson (1992). The PNS 
also spins up as a result of angular momentum conservation. 
The hydrodynamical power of the wind decreases with the neu- 
trino luminosity and there is a point in time when the magnetic 
energy density begins to dominate the thermal and kinetic wind 
energy density. An Alfven point will form, and efficient ex- 
traction of the rotational energy from the PNS will commence. 
The mass loss rate wiU Ukely increase in this phase despite the 
decrease in the neutrino luminosity as Ra approaches ^^sonic, as- 
suming that Ri,il exceeds Cs(R,y) (see §5). This is a result of 
the exponential dependence of M on fi^ in the centrifugal hmit. 
For a rapidly rotating PNS, our calculations show that for rj 
close to 2, of order 10^^ erg is naturally extracted in the first 
few seconds. In §5.3 we showed that even in the case ry = 3, 
by appealing to slightly higher Bq and fi, of order 10^^ erg 
can be extracted on a few-second timescale. The velocity of 
this wind material evolves quickly to near c and the NMC wind 
drives a secondary shock into the post-supernova-shock mate- 
rial. This secondary, but much more energetic shock encounters 
the slower supernova shock while still inside the progenitor of 
a Type-II or Type-Ibc supernova (probably at a radius < 10^ 
km).'^ In this way, a hyper-energetic supernova is created. 



Only in a very compact progenitor will the secondary wind 
shock collide with the preceding supernova shock outside the 
star. This could plausibly occur in accretion-induced col- 
lapse of a white dwarf. Otherwise, we expect this very en- 
ergetic shock to alter the nucleosynthetic signature of the hy- 
pemova with respect to normal supernovae, producing excess 
nickel as inferred from, e.g. SN1998bw and SN2003dh (Naka- 
mura et al. 2001a,b; Maeda et al. 2003; Woosley & Heger 
2003). If magnetars are bom in progenitors with extended en- 
velopes, these strong magnetocentrifugal outflows might also 
prevent late-time fallback accretion (Chevalier 1989; Woosley 

However, for a given increase in Bq, rf decreases more than rf^. We find that 

Duncan & Thompson (1992) showed that magnetar spindown via vacuum dipole radiation can power hyper-energetic supernovae in Type-II progenitors, but the 
spindown timescales via this mechanism are not fast enough to explain hyper-energetic Type-Ibc supernovae. 



cm s~ . Thus, even in the more pessimistic dipole limit (77 = 3) 
a significant amount of rotational energy can be extracted on 
a timescale much shorter than that inferred from a naive and, 
in this very early phase inappropriate, application of vacuum 
dipole spindown. 

Interestingly, if Bq is increased, Ra increases, thus decreasing rf^ 
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& Weaver 1995; Fryer et al. 1996). Although the simple anal- 
ysis presented here assumes sphericity, because the outflow 
is magnetocentrifugally driven, we expect asymmetries in the 
ejecta due to colUmation (e.g. Sakurai 1985; Begelman & Li 
1994). 

6.2. Gamma Ray Bursts 

Although the time evolution is complex and depends sensi- 
tively on Bo, our models generally achieve large F as M de- 
creases at /?L (see middle lower panel in Fig. 3). Even though 
the velocity of the matter at /?l is only mildly relativistic, 
asymptotically the energy of the field may be transferred to the 
wind material, yielding ultra-relativistic velocities. If all of the 
magnetic energy at /?l is eventually converted to kinetic en- 
ergy, the asymptotic limiting Lorentz factor (700) of the matter 
is r, given by eq. (20)."' Our exploration of the spindown of 
PNSs with 1 ms initial periods shows that F may be very large 
with E in the range > 10^' erg s"^ In fact, we obtain of or- 
der ~ 10^' - 10^^ erg emitted on a 10- 100 second timescale 
with r > 100 (Fig. 3). This conclusion is particularly robust 
for the optimistic monopole-like field geometry with 77 = 2. The 
correspondence between the energy budget, energy loss rates, 
and Lorentz factors inferred from observations of GRBs and 
the models presented here support the conclusion that a class of 
PNSs, born highly magnetic and rapidly rotating within com- 
pact progenitors, may be responsible for cosmological gamma 
ray bursts (Usov 1992; Thompson 1994). 

The now-robust association of long-duration GRBs with 
highly energetic supemovae, as evidenced by SN1998bw 
(Galama et al. 1998) and SN2003dh (Stanek et al. 2003; Hjorth 
et al. 2003), follows naturally from our discussion of hyper- 
novae in §6.1. As described in §6.1, we can expect complex 
shock structure evolution as the energy loss rate changes and the 
outflow becomes increasingly more relativistic. The dynamics 
of the supernova-shock/wind-shock interaction with the enve- 
lope and the associated instabilities might be important both in 
producing internal shocks and time variability and ultra-high- 
energy cosmic rays (Arons 2003). There are other ways to pro- 
duce time variability. We find that large changes in F naturally 
occur as evolves with fl (see Fig. 3). This order unity vari- 
ability in F provides a mechanism for producing internal shocks 
(Thompson 1994; Rees & Meszaros 1994). Shocks in the flow 
can also be created by fluctuations in mass loading caused by 
sudden changes in M. Closed loops that trap wind material 
on the surface of the PNS may be sheared open on millisec- 
ond timescales due to convection or differential rotation or may 
be opened by continued neutrino heating of the trapped mat- 
ter (Thompson 2003). Finally, relativistic MHD instabilities 
(e.g. reconnection) may produce time-variability in Poynting- 
flux dominated outflows. 

The duration of gamma ray bursts in the magnetar model 
presented here is somewhat uncertain and could be set by ei- 
ther the spindown time of the PNS (tj) or the PNS cooling 
time tkh-'^ The latter would be appropriate if tj > tkh and 
if the magnetic field significantly decreases after the cooling 
epoch, leading to a decrease in wind luminosity and thus an 

Analyses of ideal relativistic winds predict that 700 ~ r'/' (Michel 1969; Goldreich & Julian 1970; Besldn et al. 1998). However, modeling of the equatorial 
relativistic wind of the Crab pulsar (Kennel & Coroniti 1984a,b; Spitkovsky & Arons 2003) suggests that 7cx, ~ T and not that 700 ~ T'/^, in contradiction with 
theoretical expectations. 

In principle the duration could also be set by much more 'messy' physics such as the detailed emission model and how it depends on E and F in the wind. 



end to a detectable GRB. The field at the Alfven point (rel- 
evant for the spindown power) could decrease either because 
the surface field decreases in strength as convection in the PNS 
ceases (Thompson 1994) or because the field may transition 
from monopole-like to dipole-like as M decreases. 

A further important feature of the magnetar model presented 
here is that in the monopole case, E decreases exponentially 
as a function of t with timescale given by eq. (23). This is a 
prediction of the monopole-hke magnetar model and it stands 
in sharp contrast to the energy loss rate predicted from both 
coUapsar-type models (Woosley 1993) whose E must evolve 
on the timescale for accretion onto the central compact object. 
Equation (23) also shows that if the field is increased to lO'* 
G that r]^ ~ 0.3 seconds, which is in the middle of the distri- 
bution for short-duration GRBs. Similarly, allowing for lower 
Bq and f2, we further speculate that if magnetars are responsi- 
ble for long-duration GRBs, tiiat much longer (> 1000 s) GRBs 
may exist and could be detected by upcoming observations with 
improved sensitivity (e.g. SWIFT). If one allows for a dynamo 
mechanism operating efficiently within the convective core of 
protoneutron stars (Duncan & Thompson 1992; Thompson & 
Duncan 1993) then Bq is directly proportional to Q to some 
power, depending on the details of the dynamo. If there is such 
a relationship between Bq and fi, then the diversity of GRB 
sources, both long and short GRBs (KouveUotou et al. 1993) as 
well as hyper-energetic supemovae, may be grouped in a one- 
parameter family, a function only of rotation period. In this 
way, the birth of rapidly rotating PNSs may be responsible for 
a diverse variety of astrophysical explosions. We speculate that 
short bursts come from more rapidly rotating progenitors with 
no extended envelope. For example, the short GRB popula- 
tion may be due to accretion-induced collapse of a white dwarf 
(leading to a sub-millisecond PNS), whereas the long popu- 
lation is from collapse of Type-Ibc progenitors (perhaps only 
~ 1 ms periods). Another possibiUty for short GRBs is that the 
merger of a neutron star binary creates a black hole or massive 
rapidly rotating neutron star that drives a wind similar to those 
detailed in this paper (Thompson 1994). The scaling relations 
derived here are appUcable in this context as well for magneto- 
centrifugally driven flows off the accretion torus produced from 
the inner regions of a black hole accretion disk. In either case, 
the short bursts would trace an old stellar population, whereas 
the long bursts would be associated with star-forming regions. 

That the birth of neutron stars with magnetar-like field 
strengths might be the astrophysical origin of GRBs is not 
new. The idea presented here follows that of Usov (1992) and 
Thompson (1994) and touches on any model of GRBs requir- 
ing Poynting-flux dominated outflows (e.g. Lyutikov & Bland- 
ford 2003). In Usov (1992) and Thompson (1994), only a pure 
dipole field geometry was considered. In this paper we find po- 
tentially much larger energy loss rates for a given Bq and O as 
a consequence of different field geometries. In particular, mo- 
tivated by the literature on stellar winds and their angular mo- 
mentum evolution, we have considered monopole field topolo- 
gies in the spirit of Weber & Davis (1967). Further considera- 
tion of rapid rotation coupled with high mass loss rates (see §3; 
Mestel & Spruit 1987) supports this view. The energy extrac- 
tion rates and braking suggested by our pure monopole calcu- 
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lations (see §4 and §5.3) are likely overly optimistic. However, 
we await multi-dimensional MHD models of relativistic out- 
flows. 

In addition to the issue of field topology there are other im- 
portant potential difficulties inherent to this GRB mechanism. 
Similar to the collapsar model of GRBs (Woosley 1993; Mac- 
fadyen & Woosley 1999), the energy injection rate must remain 
high on the timescale of the relativistic ejecta to break out of the 
progenitor. Given the timescales we have derived and tkh, this 
is plausible for compact Type-lbc core-collapse scenarios. In 
addition, we rely on conversion of magnetic energy into kinetic 
energy asymptotically, a process that is not well understood in 
models of pulsar winds. Issues of colHmation and - possibly - 
jet formation require multi-dimensional models and we assume 
that many of the structures generic to collapsar jets will obtain 
here as weU (Aloy et al. 2000; Zhang et al. 2003), altiiough flie 
flow wiU be Poynting-flux dominated. 

6.3. Conclusions 

We have estimated spindown timescales for protoneutron 
stars born rapidly rotating and with large surface magnetic 
field strengths. The combination of magnetocentrifugal ef- 
fects and neuttino energy deposition combine to produce a 
wind capable of significantly effecting the early angular mo- 
mentum evolution of newly born protoneutron stars. In ad- 
dition, we have constructed steady-state models of neutrino- 
magnetocentrifugal winds, including magnetic fields impUcitiy 
and approximately by forcing the wind matter into sttict coro- 
tation with the PNS. 

We find that magnetar-like surface fields (~ 10^^ G) dom- 
inate the thermal pressure and kinetic energy density of PNS 
winds in the first few seconds after the preceding supernova for 
a wide range of initial spin periods. We show that as the mass 
loss from the PNS abates, the outflow becomes increasingly 
faster, transitioning from non-relativistic to relativistic outflow 
in the first few to tens of seconds. We explore this progression 
and provide evolutionary models that bridge the transition. We 
find that non-relativistic spindown timescales can be consid- 
erably shorter than those inferred from simple vacuum dipole 
spindown. This results both from physics of the Alfven point 
and the large enhancement in mass loss rate expected in models 
with spin periods less than '--^ 2 ms (see §5 and Fig. 2). Moti- 
vated by the Uterature on stellar winds from magnetic rotating 
stars, we have carried out calculations for monopole-like field 
geometries in addition to pure dipole fields (see §3 and §5). 
Although optimistic in computing spindown, quasi-monopolar 
fields are physically motivated by the large mass loss and rapid 
rotation of the systems considered. We are currently working 
on axi-symmetric two-dimensional MHD models to more fully 
explore the issue of global magnetic field topology and angular 
momentum loss. 



Regardless of field topology, we find generically in our mod- 
els that if magnetars are born with 1 ms initial periods and 
Bo ~ 10'^ G that more than 10^' erg of rotational energy can 
be extracted in the first few seconds and that the velocity of 
this matter will approach c as the Alfven radius approaches the 
light cylinder. As the mass loss rate continues to plummet, the 
asymptotic limiting Lorentz factor of the ejected matter may 
approach several hundred. We therefore propose these objects 
as candidate central engines for hyper-energetic supernovae and 
cosmological gamma ray bursts of both the long and, possibly, 
the short variety (see also Usov 1992, Thompson 1994, and 
Lyutikov & Blandford 2003). 

Although more mundane from an energetic point of view, if 
monopolar field geometries are at all relevant for slow rotators 
(P > 20 ms), efficient angular momentum loss results. We find 
that if magnetars are born with the spin periods of normal neu- 
tron stars, taking ry = 2, that a PNS can be spun down from a 
spin period of ^ 50 ms to a spin period of ~ 1 second in ~ 10 
seconds of wind evolution. 

We emphasize that for any field strength there is a point in 
time during the cooling of the nascent neutron star when the 
magnetic field dominates the wind dynamics. If the field is 
strong enough, this happens early in the wind phase and effi- 
cient extraction of rotational energy is likely a natural conse- 
quence. The neutrinos still carry away the neutron star's bind- 
ing energy (~ 3 x 10^^ erg), but their energetic coupUng with 
the surrounding environment is (necessarily) weak. In marked 
contrast, the rotational energy (although only of order 10^^ erg 
for a fast rotator) may be efficiently coupled to the wind ma- 
terial via the strong magnetic field. It is in this way that the 
wind's energetic presence within supernova remnants and in the 
spin distribution of neutron stars may be observed. 
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Fig. 1. — Thermal pressure (P, solid line), kinetic energy density (pv^/2, short dashed Une; long dashed line), and magnetic energy density (B^/Stt, dotted 

lines) for Bq = l.6x lO'^ G and Bq = 7.5 X lO'^ G for r; = 2 and rj = 3, respectively, as a function of radius. The PNS has spin period of 1 ms (SI ~ 6283 rad s"', 
Rl — 47.7 km is the corresponding Ught cylinder radius), Lp^ = 4 X 10^' erg, Ri, = 10 km, andM= 1.4 Mq. The mass loss rate isM ~ 1.6 x 10"' Mq s"' (compare 
with Fig. 2). Model numbers were chosen to closely follow those employed for the scaUng relations in §4. 



13 



T 1 1 1 1 — I I I I 1 1 1 1 1 — l—TT 




100 1000 104 

n (rad s-0 



Fig. 2. — Log of the mass loss rate as a function of f2 for L^^ = 8 x 10^' erg s ^ and Lp^ = 0.5 x 10^' erg s ^ assuming corotation with the PNS surface, = rO. 
(soUd Unes). Also shown are results assuming no corotation (dashed Unes), but taking = RvCl(Rv/r), appropriate for angular momentum conservation. For 
R^Cl > Cj(ff^), M increases approximately exponentially with Cfi in the models assuming corotation. 



14 




t (s) I (s) t (s) 



Fig. 3. — Representative time evolution for the spin evolution of a rapidly rotating highly magnetic PNS with a monopole field geometry (tj = 2). Calculations 
with four different surface magnetic field strengths are shown: Bq = 3 X 10'^ G (A), 1.5 X 10'^ G (B), lO'^ G (C), and 7 X lO''* G (D). The upper and lower right 
panels show the mass loss rate (M [Mq s"']) and the total neutrino luminosity (logjgLi/ [10^' erg s"']), respectively. The spin period (P [ms]; upper left panel), 
limiting asymptotic Lorentz factor (r; lower middle panel), spin energy extracted (E [\CP^ erg]; lower left panel, solid lines), and energy extraction rate (E [lO'' 
erg s"']; lower left panel, dashed lines), are shown. In addition, the upper-middle panel shows i?L (solid lines) and Ra as computed with ri = 2 (dashed lines). For 
comparison, in the plot of F versus t, we also include the evolution for Bq = 1.3 x 10^' G (dotted line), which falls between models (B) and (C). Note the interesting 
time evolution of F as M drops exponentially after ? = 40 s, in response to the assumed change in L„ . 
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Fig. 4. — Representative time evolution for the spin evolution of a slowly rotating highly magnetic PNS with a dipole field geometry (rj = 3). Spin period (P [ms]; 
upper left panel), limiting asymptotic Lorentz factor (F; lower right panel), Ught cyhnder radius (Sl [km], upper right panel, soUd line), Alfven radius (R/^ [km], 
upper right panel, dot-dashed hne), spin energy extracted {E [10^' erg]; lower left panel, solid hne), and energy extraction rate (£ [10^' erg s"']; lower left panel, 
dashed line), are shown. As an extreme example, Bq was taken to be 5 x 10'* G. 
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Fig. 5. — Representative time evolution for the spin evolution of a rapidly rotating liiglily magnetic PNS witfi a dipole field geometry (t] = 3). Calculations with 
three different surface magnetic field strengths are shown: Bq = lO"" G (A), 5 X lO" G (B), and lO'^ G (C). Panels, units, and linestyles are the same as in Fig. 4. 



